#!/usr/bin/env python

from math import e, sin, cos

def exact_solution(x):
	c1 = -0.5
	c2 = 1./(2.*sin(20.)) * (cos(20.)- 1./e**20.*(1.+20.)**2 )
	return (x*x-1.)*sin(x)
#	return c1*cos(x)+c2*sin(x)+0.5*1./(e**x)*(1.+x)**2

